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ABSTRACT 

Dust production among post-main-sequence stars is investigated in the Galactic globular cluster 47 
Tucanae (NGC 104) based on infrared photometry and spectroscopy. We identify metallic iron grains 
as the probable dominant opacity source in these winds. Typical evolutionary timescales of AGB 
stars suggest the mass-loss rates we report are too high. We suggest that this is because the iron 
grains are small or elongated and/or that iron condenses more efficiently than at solar metallicity. 
Comparison to other works suggests metallic iron is observed to be more prevalent towards lower 
metallicities. The reasons for this are explored, but remain unclear. Meanwhile, the luminosity at 
which dusty mass loss begins is largely invariant with metallicity, but its presence correlates strongly 
with long-period variability. This suggests that the winds of low-mass stars have a significant driver 
that is not radiation pressure, but may be acoustic driving by pulsations. 

Subject headings: stars: mass-loss — circumstellar matter — infrared: stars — stars: winds, outflows 
— globular clusters: individual (NGC 104) — stars: AGB and post-AGB 



1. INTRODUCTION 

Stellar mass loss is of critical importance to the later 
stages of stellar evolution. Stars evolving up the red giant 
branch (RGB) lose mass via stellar winds. The remain- 
ing mass of a star leaving the RGB tip is the primary de- 
ter minant of its position on th e horizontal branch (HB; 
e.g. [Roodl [1971 ICatelanI [2000l) . Subsequent mass loss 
on the asymptotic giant branch (AGB) ejects the star's 
entire hydrogen envelope, creating a post-AGB star and 
(perhaps) a planetary nebula (PN). The timescale and 
end-point of AGB evolution may therefore not be de- 
termined by the nuclear burning rate, bu t by mass loss 
(|van Loon et al.lll999t [Bover et al.ll2009aD . Winds from 
these stars are the most significant producers of inter- 
stellar dust and therefore provided the heavy elements 
present in Population I stars, including t he Sun and So- 
lar System (|Gehrzlll989l: ISedlmavdfT99l IZinneri[200l . 

Debate exists around many finer points of the mass 
loss process, particularly its variation with fundamental 
parameters such as luminosity, temperature and metal- 
licity. Firstly, the amount of mass loss that occurs 
on the RGB remains largely unmeasured. RGB mass 
loss appears to determine the 'second parameter' (af- 
ter metallicity) that defin es the morphology of HBs in 
globular clusters (GCs'l (iFusi Pe cci & Bella zzinil 119971: 
iLamers fc Cassinellil [19991 : ICatelanI I200a) . Secondly, it 
is not proven where dust formation stars, and whether 
RGB stars can form significant amounts of dust in their 
ejecta. If many RGB stars do produce dust, it would 
add a large number of dust factories ejecting dust into 
a galaxy's interstellar medium (IS M) (cf. the missin g 
dust source in many galaxies; e.g. Matsuura et al.|[2009[ ). 
lOriglia et al.l (|2007| ) suggest that dust production oc- 
curs episodically over a large portion of the RGB, while 
[McDonald et al.l (|2009( ) have presented evidence that 
dust production is confined to the RGB tip, though 
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gaseous m ass loss is still present all along both the RGB 
and AGB ([McDonald fc van LoonI 120071 : iMeszaros et all 
[200I IDupree et al.ll2009D . Dust can also help d rive a 
wind through radiation pressure on grains (e.g. [Lewis! 
Il989f l. though it is not certain whether winds which lack 
high-opacity amorphous carbon can be accelerated in 
this manner (IWillso nI ^M)?! IWoitkd [2nM IHofned [MF^ 
iLagadec fc Ziilstral 12008: M cDonald et~al|[2010[) . Dust- 
less stars typically lose mass via magneto- acoustic pro- 
cesses in their chro mospheres (iHa rtmann fc MacGrcgo3 
[1980; 'Dup ree et an[l98l iMcDwiald fc van LoonI I200C 
Finally, the composition of circumstellar dust and the 
chemical processes by which it forms are only partly de- 
fined. Of particular interest is the role that metallic iron 
grains play in dust formation, and ho w metallicity af- 
fects dust production and composition ([McDonald et al.l 

[2oTnh . 

In this study, we aim to address these points. To do 
this, we explore dust production in the Galactic glob- 
ular cluster 47 Tucanae (NGC 104). Globular clusters 
are unique laboratories with which we can study late- 
stage evolution of low-mass stars: globular clusters con- 
tain stars that are among the oldest in the Universe, 
and most clusters are comprised of a single (or at least 
dominating) stellar population with little internal vari- 
ance in age or chemical composition. Comparisons be- 
tween clusters therefore probe differences in mass loss 
created by metallicity, while comparisons within clusters 
probe differences created by stellar temperature and lu- 
minosit y. Throughout th i s pape r, we will refer to similar 
studies: [McDonald et al.l ( |2009l ). which cove rs the cluster 
wCentauri (NGC 5139; hereafter Paper V): iBover et al.l 
(|2009H), which covers the cluster NGC 362 (hereafter Pa- 
per II); and a work which accompanies this paper (sub- 
mitted ApJS; hereafter Paper III). 

The cluster 47 Tuc itself is among the Galaxy's 
most massive and, co-incidentally, one of the clos- 
est to us. Its distance h as been est i mated photo- 
metrically as 4190-4700 pc ([Harrisl [l99l iPercival et al.l 
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McDonald et al. 



I2002t iSalaris et al 



(jMcLaughlin et al 



mo 



200' 



and kinematically as 4000 pc 
with 4500 pc being the median 
distance determination in the literat ure. Its mass is 6- 
9 X 10^ (ISco tt & Rose 1975; Mc vlan fc M"avoi1[l985l : 
IMandushev et al.lll99lO . Its metallicity is determined to 
be roughly [Fe/H] = -0.7, though various studi es have 
placed it anywhere from fFe/H| = -0.89 to -0.66 ([Harris 

i[ 



1996HCarretta fc Gratton|ll997i:lMc"William fc Bernstein 

m 



2008t iWorlev et al.l 120101 ) . The cluster suffers from 
littl e inter s tellar reddening {E{B — V) = 0.04 mag 

— [Harris! I1996D and is estimated to be oy er 11 

Gyr in age (iC habovcr et al.' '1992*: 'Zoc cah et all 120011: 
Grundahl et al.ii2002^,Percival et al.., 20021: iGratton et al.l 



20031: IKaluznv et all 120071 ISalaris et"aT] 120071) . The 
cluster hosts an ionised intra-cluster medium (ICM; 
iFreire et al.ll200ll ). This ionised medium is undetected 
in any other cluster (M15 has a neutral medium with 
a dusty component: lEvans et al.l 120031 Ivan Loon et al.l 
l2006bl : iBover et al.l I2006D . Integrating the observed in- 
tracluster electron density over the cluster core yields 
~0.1 Mq of H II, represent i ng >5% of the in tracluster 
hvdrogen (jSmith et al.l[T990l : iFreire et al.ll2001| ). The to- 
tal stellar mass loss from 47 Tuc's stars therefore gives 
us insight into the fate of gas and dust expelled by stars 
in this cluster. 

Paper HI more-precisely examines which stars have 
mid-infrared excess, indicative of dust, finding again that 
we cannot corroborate the presence of dust around the 
fainter giants. In this work, we use the data and results 
from Paper HI to derive mass-loss rates for individual 
stars and (by summation) the cluster as a whole, and in- 
vestigate the composition of the circumstellar dust pro- 
duced by analysing archival mid-infrared (mid-IR) spec- 
tra. 

2. MASS LOSS AND DUST PRODUCTION 

2.1. The stellar mass-loss rate in J^l Tuc 

2.1.1. The input data 

The data for this study come from the accompanying 
Paper HI, which identifies those stars in 47 Tuc which 
have infrared excess, indicative of circumstellar dust. 
This list of 22 objects is an exhaustive list of those stars in 
the cluster which we can be relatively certain are produc- 
ing dust, given currently available data. These objects 
are listed in Table [TJ which also includes the parameters 
we determine in later sections. 



2.1.2. A measure of mass-loss rate 

Accurate — or even approximate — mass-loss rates 
are notoriously difficult to measure. One must assume 
that the wind has certain properties. Perhaps the most 
crucial of these is the relationship between wind veloc- 
ity and distance from the star. We here assume a con- 
stant, spherically-symmetric, purely radiatively-driven 
wind, starting from a near-statio nary, dust-forming layer . 
As discussed previously (Paper I: iMcDonald et al.ll2010l ). 
this may not be the most appropriate description, but 
suffices due to lack of contrary evidence: we return to 
this point in the discussion. 

^ Th e updated Harris cluster catalogue can be found at: 
|http: / /www.physics.mcm aster.ca/~harris/mwgc.dat 



Before calculating an estimated mass-loss rate, we can 
create a measure of the amount of circumstellar dust sur- 
rounding a star, irrespective of most modelling uncertain- 
ties. To do this, we define a quantity /r, the fraction of 
the star's luminosity that is reprocessed by circumstel- 
lar dust. For low optical depths, this should be directly 
proportional to the surface area of dust grains, therefore 
approximately proportional to the mass of circumstellar 
dust. The value of /r can be determined by taking the 
difference between the observed SED and model photo- 
spheric SED in the infrared (yielding the infrared excess) 
and dividing this by the total photospheric luminosity 
(the value given in Paper HI). We list this quantity in 
Table [T] Errors in the value of /r are introduced by dif- 
ficulties in determining the integrated luminosity from 
photometry alone and from the uncertainty in the pho- 
tospheric spectrum (both temperature and luminosity). 
The former can be approximated as A/r = y/l -\- /r — 1, 
or roughly 50% of /r, but is dominated by photometric 
uncertainty at short (2.2-10-/im) wavelengths. The lat- 
ter gives a constant uncertainty of around A/i- w 0.4%. 
We reduce this latter error by fixing the photospheric lu- 
minosity to the if -band flux for all stars except V26 and 
x03, which are fixed to the 3.6-/j,m flux. The value /r 
could not reliably be measured for several stars with low 
mass-loss rates. 

The stars appear to fall into three main categories: 
those with large (~5%) reprocessed flux (Vl-4, V8), 
those with low (^2%) reprocessed flux concentrated to- 
ward long wavelengths (V13 and VIS), and those stars 
with low (<2%) reprocessed flux concentrated toward 
shorter wavelengths (the remainder). 

2.1.3. Defining the mass-loss rate 

Given the limitations of our understanding and of the 
model presented above, we stress that the mass-loss rates 
and velocities we derive are based only on one model 
for wind acceleration, but should be comparative both 
internally and with other works which use these same 
assumptions. 

Under these assumptions, we calculate mass- loss rates 
by modelling the SED using the radiative transfer pack- 
age, DUSTY (Nenkova et al. 1999). As inputs, we use our 
interpolated MARCS model of the stellar s pectrum, and 
a sta ndard MRN grain size distribution (jMathis et al.l 
Il977f) . Within our DUSTY model, we vary the dust com- 
position, temperature in the dust-forming region (the in- 
ner edge of the dust envelope) and optical depth of the 
wind. We change these parameters iteratively until a 
best fit is achieved. 

We relate the dusty mass-loss rate (Mdusty) to the 
actual mass- loss rate (M) as pr eviously (Paper I; Paper 
II), using the following formula ()van Loonl[2000l) : 



M = Mdusty 



L 



,104 Lq 
lO-^AfouSTYi/Wo 



3/4 



200 



P 



ifi 3 g cm 



1/2 



[Mq yr- 



(1) 
(2) 



for luminosity L, dust-to-gas ratio ^, grain density p, and 
terminal wind velocity v^o ■ We assume that p = 3g cm~^ 
for silicates, 4 g cm~^ for aluminium oxides, 5.7 g cm"'^ 
for FeO, and 7 g cm~^ for metallic iron. For these three 
species, we use the optical constants from iDraine fc Lei 
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TABLE 1 

Calculated wind properties and literature pulsation properties for potentially dusty stars in 47 Tuc. 
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(mag) (km /s) (days) 






Ks 


[3.6] 


[8] 


(/oj (lU 
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Sil. Fe 


AI2O3 


FeO 




V 1 


4.03 


20 


221 


3623 


4824 


0.936 


1.173 


2.213 


3.6 


21.0 


4.0 


900 


12 88 










V8 


1.6 


16 


155 


3578 


3583 


0.904 


1.210 


2.135 


5.9 


14.7 


4.0 


900 


7 85 


7 


1 




V2 


2.78 


23 


203 


3738 


3031 


1.103 






4.3 


12.0 


3.8 


900 


5 95 








5 


V3 


4.15 


22 


192 


3153 


2975 


0.954 


1.198 




3.2 


9.4 


3.2 


1000 


100 








6 


V4 


1.5 


18 


165 


3521 


2603 


0.932 


1.167 


1.902 


5.6 


11.7 


3.6 


900 


2 95 


2 


1 


~RGB-tip 


V26 


0.3 






3500 


2541 


1.073 


0.807 


0.995 


0.9 


5.9 


3.8 


1000 


100 






7 


LWIO 


1.2 




221; 


3543 


2324 


1.073 


0.876 


1.017 


0.6 


4.2 


3.2 


1000 


100 










V21 


1.0 


7 


76: 


3575 


2301 


0.977 


1.083 


1.355 


1.1 


4.9 


3.6 


1000 


93 


5 


2 




LW9 


1.1 




74 


3374 


2204 


0.973 


1.003 


1.228 


1.1 


6.0 


3.3 


950 


100 








V27 


0.3 




69 


3374 


2140 


1.023 


0.929 


1.218 


1.4 


3.3 


3.6 


1050 


100 








Leel424 ■■■ 




65 


3565 


2122 


1.006 


1.182 


1.390 


>4 


4.4 


4.0 


1050 


100 








A19 


0.7 




60 


3526 


2096 


1.034 


0.992 


1.240 


0.8 


3.6 


4.2 


1100 


100 








LW12 


0.68 




116 


O 1 lO 


9070 


1.082 


0.967 


1.093 


\ -L 
















x03 








3816 


1640 




0.958 


1.156 


3.2 


5.1 


1.8 


750 


100 






11 


LW19 


0.1 




40 


3738 


1638 


1.078 


0.991 


1.080 


<1 


low 














V20 


0.8 




232: 


3602 


1575 


1.053 


0.928 


1.039 


<1 


low? 












10,11,12 


V22 


0.5 




62 


3684 


1528 




2.504 


2.834 


<1 


low 












7 


V6 


0.6 


7 


48 


3763 


1374 


1.001 


1.000 


1.167 


<1 


low 














V5 


0.7 


8 


50 


3741 


1363 


1.003 


1.100 




<1 


low 
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V18 


0.3 


5 


83: 


3692 


1297 


0.989 


1.004 




2.4 


5.3 


2.4 


800 


33 67 








11 


V23 


0.5 




52 


3775 


1259 


1.023 


0.985 




<1 


low 












7 


V13 




12 


40 


3657 


1029 


0.990 


1.106 




1.2 


4.1 


2.2 


800 


85 


15 


5 
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^Optical and radial velocity variability (optica l in V-ba n d, wh er e known) w it h per io ds (colons denote ap proximate v alues, preference 
given to the longest specified period), from IClemen^ II 19971 ): IPoimanskil 120021) : ILebzelter &c Wood! 12005?) : and ILebzelter et al.l 
(2005). ^Total (gas + dust) mass-loss rate and terminal velocity assuming a purely radiatively-driven wind, assuming ^ = 1/1000. 
^Temperature of the inner edge of the circumstellar dust envelope. "^Implied fraction of grains in 'astronomical' silicates, metallic iron, 
aluminium oxide and iron oxide. ^No Spitzer IRAC data, fit based on literature photometry and MIPS 24-^m data. ^Unusually 
cool object. ^ Photometry uncertain du e to blending, fits made to Spitzer IRAC data only. ^May switch between two states of 
variability, see ILebzelter fc WoodI 1)20051 '). ^Values titled low are undetermined, but liable to be < 2 X 10~^ Mq yr~^. '^''8- & 24-/im 
emission only, suggesting the presence of silicate grains. ^^Wind may be impossible to sustain (see text). ^^Very poor-quality JHKs 



photometry, IR excess difficult to determine. 



(|1984D : lBegemann et al.l (|1997| ): and lOrdal (|1988f ). 
respectively. The errors in the above assumptions are 
substantial, and are discussed in Section [2.1.5l 

The dust-to-gas ratio, i/i, can be estimated as: 



1 
200 



X lOlF^/Hl 



(3) 



which is approximately 1/1000 for 47 Tuc. This assumes 
that the fraction of condensed metals is identical for all 
metallicities. We discuss the accuracy of this assumption 
in Section [^.1.5 1 

2.1.4. The stellar mass-loss rate 

Table [T] shows our calculated mass-loss rates and in- 
ferred wind parameters and dust compositions along- 
side literature pulsation properties. The striking fact 
about the dust composition is the fraction of iron grains 
present. The presence of metallic iron in these stars is de- 
batable, as its featureless infrared emission can be mim- 
icked by many other phenomena, including large silicate 
grains and emission from thi ck circumstellar molecular 
layers (jMcDonald et al.ll2010[ ). We consider iron to be 
the most likely source of emission, though the fractions 
of metallic iron listed in Table [T] may not be particularly 
accurate, and we discuss this further in §3.51 

Many of the mass-loss rates are too low to be sensi- 
bly measured: we can merely say that some circumstel- 
lar dust probably exists. Those rates that can be mea- 
sured have their SEDs shown in Figure [TJ The photom- 



etry is taken from the sources cited in Paper IIF|- Mid- 
IR spectra are includ ed from Ivan Loon et al.l (|2006aD 
(hereafter vLMO), and ILebzelter et al.l (|2006f ) (hereafter 
LPH). These objects are summarised in the following 
luminosity-ordered list: 

• VI: This star is clearly identified as the clust er's dom- 
i nant dust producer by vLMO, LPH and llta et all 
()2007f ) (hereafter ITM) . VI has clear excess at all wave- 
lengths > 3^m. This excess, which we attribute to 
metallic iron dust, dominates even over its strong sili- 
cate feature. 

• V8: No silicate feature was seen in this star by vLMO, 
but a clear silicate and aluminium oxide feature was 
seen by LPH. It has previously been suggested (LPH; 
Paper I) that V8 may have variable mass loss. How- 
ever, the AKARI photometry at 3, 4, 7, 11, 18 and 23 
/Ltm shows no variation in the (dominant) iron contri- 
bution. This means that if the mass loss is variable, 
strong variation only occurs in the silicate production. 
V8 has a complex IR spectrum showing peaks from 
silicates, aluminium oxide, and from crystalline min- 
erals at 9.7, 11.3, 13.1, 20 and 32 ^ m. We discuss 
the carrier (s) of these bands in ij3.5l Note that the 
amplitude of the 10-/j,m silicate peak is substantially 



_f Optical data: IStetsonI l|200a 'l: IZaritskv et al.l 

alaris ct al. ( 2007|) ; near-IR data: Skruts kie et al.l 
Salaris et al. (200^; mi d-IR data : [Ori gl ia ct al.l 11997] , 
van Loon et al. (2(W6i); llta et aU 120071): iBarmbv et aU 
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Gordon et al. (2009); Hover et al.l 1)2013 ): llshihara et al.l HgOlQ ) 
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Fig. 1. — IR portions of the SEDs of dusty stars, plotted in Ra yleigh -Jeans units. Large (red) and small (black) points show literature 
photometric and smoothed spectroscopic data, respectively (see ^2.1.41 for reference list). Spectroscopic data have been scaled in flux to 
match the photometric data. The long-dashed (green) and short-dashed (blue) lines show? the modelled photospheric contribution and total 
flux, respectively. 



over-estimated in our model compared to the 20-/im 
silicate peak. This implies a lower temperature for the 
silicate grains than the metallic iron grains, which we 
also discuss in 



V2: A silicate dust feature was observed by LPH, 
though it was too weak to have been seen by vLMO. 
Iron dominates the IR flux, however. The 10- and 20- 
/im silicate peaks are not well modelled, as in V8 (see 
also 3331). 



• V3: No silicate emission was seen by either LPH or 
vLMO, but strong, featureless (metallic iron) emission 
is present, which is confirmed by AKARI photometry. 

• V4: The spectrum in LPH shows complex dust miner- 
alogy, akin to that in V8. This star is roughly at the 
RGB-tip, as estimated from the luminosity function. 

• V26: Blended with a star half its brightness, 2.2" to 
the east, the near-IR, AKARI and 24-/im photometry 
are affected. The amount of IR excess is difficult to 
determine here, but is likely to be quite small. No 
discernible features are present in the spectrum from 
vLMO. 

• LWIO: Also within the dense cluster core, it appears 
that the AKARI photometry (at 9 & 11 nm) for this 



object is affected. No features are detectable in the 
spectrum from vLMO. 

• V21: Observed by LPH, the spectrum shows the same 
features as V8. These features, and the silicate contri- 
bution, are considerably weaker in V21. 

• VI 1: No obvious features were seen in the spectrum 
of vLMO, this star was only covered by Spitzer at 3.6 
and 5.8 /im. AKARI 9-/im photometry suggests it 
may have some excess there, but both the Spitzer 8- 
/im photometry and the vLMO spectrum (which covers 
9-/im) show no excess. The derived mass-loss rate can 
probably be treated as a rough upper limit. 

• LW9 & V27: An excess consistent with circumstellar 
dust is present in both the Spitzer and AKARI pho- 
tometry. The chemistry cannot be well-constrained by 
this photometry, but it appears that the flux contribu- 
tion from components other than iron is comparatively 
small. 

• Leel424: Lacks any data redwards of 9 /im, making 
constraining the mass-loss rate, dust temperature and 
dust chemistry very difficult. The values given should 
be treated as order-of-magnitude estimates only. 

• A19: As LW9 & V27, however the AKARI 7-fim data 
suggests little to no mid-IR excess in this star. 
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• V20 & V22: not shown in Figure[TJ they have no near- 
IR photometry attributed to them, meaning it is dif- 
ficult to determine the amount of IR excess present. 
Despite this, the dust temperature is relatively-well 
constrained by the presence of Spitzer 24- /im photom- 
etry. 

• x03: Located in the cluster core, this star suffers from 
the same problems as V20 & V22, but clearly displays 
a stronger excess, making an approximate determina- 
tion of the mass-loss rate and wind characteristics pos- 
sible. 

• V18: Observed by both LPH and vLMO, this star has 
a clear silicate feature. This feature was not claimed in 
vLMO due to poor signal-to-noise, but the detection 
appears to match that from LPH well. It shows sub- 
stantial mass loss compared to other stars of similar 
luminosity. The 10- and 20-/im silicate peaks are not 
well modelled, as in V8 (see also ^?>A\ . 

• V13: Highlighted in both LPH and ITM, the spec- 
trum of this star shows an excess similar to that of 
V8. The contribution of iron dust is comparatively 
small, though would appear to be non-zero. This star 
is isolated from the cluster core, but has poor-quality 
photometry due to its comparative faintness. 

2.1.5. Error m the mass-loss rate 

The accuracy of the wind parameters we find is hard to 
determine. The major sources of error are listed below: 

1. A 30% uncertainty inherent in dusty due to an un- 
specified gravitational correction. This uncertainty 
is exceeded in stars with Voo ^2.5 km s""'^, where 
the effect of stellar gravity will change the variation 
of dust particle density with radius from the star 
(the stellar mass is unspecified in dusty). 

2. Uncertainty in the underlying photospheric spec- 
trum. As the iron dust contribution is essentially 
a modified blackbody, one can fit a dust tempera- 
ture to the 3.6- to 8-/zm data and extrapolate this 
fit backwards to ifg-band, where the dust contri- 
bution is essentially zero. The photospheric con- 
tribution can then be set from the UBVIJHKs 
photometry, but the accuracy of this fit depends on 
the stellar photometry. (This would equally hold 
for other possible dust species, e.g. amorphous car- 
bon.) 

3. Uncertainty in the amount of IR excess. Apart 
from identifying how much fiux is generated by the 
photosphere, missing photometry and photometric 
errors cause uncertainty in the combined flux of the 
(photosphere -f dust) contributions. 

4. Uncertainty in the dust composition. It is not 
clear for a number of these stars (V26, LW9, V27, 
Leel424, A19, x03 and any star marked "low") 
what the dust chemistry is. Many of these stars 
have insufficient coverage in the 9-23 /im region 
to determine whether they contain any dust com- 
ponents other than metallic iron. LWIO and VI 1 
are not covered by LPH: ground-based spectra by 
vLMO show no obvious spectral features, so we 



presume a pure metallic iron wind, but the spectra 
are of low signal-to-noise. The presence of silicates 
can significantly increase the implied mass-loss rate 
due to their low opacity. Any carbonaceous dust 
present would also be hidden by the modified black- 
body of metallic iron. Carbonaceous dust typically 
has a higher opacity: if it is present, the implied 
mass-loss rate would be decreased. 

5. Uncertainty in grain size distribution and grain 
shape (assumed here to be spherical, with radius 
a). The relationship between the grains' absorption 
cross-section to radiation (cx a?) and their mass 
(oc a'^). A high grain porosity would also imply a 
decreased dust density and therefore a higher mass- 
loss rate. Particularly, differences in grain proper- 
ties between the iron and silicate grains could ex- 
plain why we model such a high proportion of iron 
grains in the wind (see §3.5|) . 

6. Uncertainty in dust temperature. The dust tem- 
perature can only be determined to within ~100 
K. As L cx R^T'^, a 15% error in temperature leads 
to an error of up to ^60% in the flux emitted by 
the circumstellar dust, depending on the implied 
error in the radius of the dust shell. This will be 
most significant in stars that are losing the least 
mass, as their dust temperatures are more difficult 
to determine. 

7. Uncertainty in the wind driving mechanism. We 
presume here that the winds are sustained by radia- 
tion pressure alone. The very low outfiow velocities 
present in at least three cases (x03, V18 and V13) 
would suggest that there is insufficient radiation 
pressure for it to be the sole driver of mass loss. Ad- 
ditional energy input, for example from pulsation, 
would increase the wind velocity and (by virtue of 
Eq. ([2])) decrease the mass-loss rate. There is a 
related uncertainty in the dust-to-gas ratio, which 
is assumed to scale with metallicity. These uncer- 
tainties cannot be quantified statistically. 



The easiest way to see how these uncertainties affect 
a star's determined mass-loss rate is to propagate them 
through the DUSTY code. The results for VI (Table [2] 
suggest that, in absolute terms, our mass-loss rates are 
little more than an educated guess. However, since a 
large number of other studies use the same assumptions, 
one can use these rates to trace differences between stars, 
providing these properties are largely invariant between 
stars. 

The parameters listed in Table [2] have different effects. 
The grain size, for example, has little effect on the mass- 
loss rate but significant effects on other wind properties. 
Conversely, the dust-to-gas ratio have significant effects 
on the mass-loss rate and velocity but little effect on the 
other wind parameters (temperature, optical depth and 
minerology). Clearly the largest effect on the mass- loss 
rate is whether the wind has an initial velocity driven 
by, for example, pulsation. Determining the radial ve- 
locity profile of the wind would reduce the error sub- 
stantially. While this has not been done accurately for 
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TABLE 2 

Error analysis of mass-loss rates for VI. 



Parameter 


5M 


SVoo 


Parameter variation used 


Other parameters affected 


Grain size^ 


±10% 


±30% 


Low: q = 2.5, r = 1 nm, R = 100 nm 
Assumed: q = 2.5, r = 5nm, R = 250 nm 
High: q = 1.5, r = 5nm, R = 250 nm 


Tinner = 750 K, 90% iron, 10% silicates 
Tinner = 850 K, 88% iron, 12% silicates 
Tinner = 850 K, 80% iron, 20% silicates 


Grain density 


-50% 


-50% 


Low: 1.63 g cm"*^ 

Assumed (solid): 6.52 g cm~'^ 




Dust formation 


±10% 


±20% 


Low: 750 K 


T-v = 0.38, 92% iron, 8% silicates 


temperature 






Assumed: 850 K 


T-v = 0.45, 88% iron, 12% siUcates 


Dust: gas ratio 


-55% 
+ 125% 


+ 125% 
-55% 


High: 950 K 
Low: 1 / 5000 
Assumed: 1 / 1000 
High: 1 / 200 


TV = 0.52, 82% iron, 18% silicates 


Velocity distribution 


~+500% 


N/A 


Assumed: — v^o = 4 km 
High: 10 km s~^ (constant) 


Tinner = 850 K 
Tinner = 900 K 


Photosphere 


±12% 


±2% 


Assumed ± 5% in luminosity 


TV decreased from 0.45 to 0.42 




±15% 


±6% 


Assumed ± 250 K in temperature 




Calculation error^ 


±30% 








Total 


-tlx 


+N/A ^ 
~-3.5 ^ 







^The distribution of g rain radii (a) is defined by a slope of a ' between sizes r and R. "^Due to an unspecified correction for stellar 
gravity within DUSTY IjNenkova et al.|[l999 '). 



such low-metallicity, dust-producing giants, inferences 
from chromospheric line modelling ( e.g. lJudge fc Stenceil 
[1991 IMcDonald fc van LoonI 120071: IVievtes et al.ll2010D 
appears to show that wind velocities do not decrease sub- 
stantially from ~10 km in the stages shortly prior to 
dust formation. We point out, however, that the highest 
mass-loss rates (and by inference the highest wind veloc- 
ities) may be ruled out on evolutionary grounds ( ^3.5p . 

We emphasise that the errors derived above are not 
formal and are only estimates. These compare relatively 
well to the adopted error of ^5 x for silicate-dominated 
winds (e.g. iSargent et al.l[2010h . where the grain prop- 
erties are somewhat better determined due to the long 
history of their study and constraints from their spectral 
features. 

In relative terms, however, we can probably assume 
that the star-to-star variations are much smaller. In par- 
ticular, the driving mechanism and hence the wind's ini- 
tial velocity should not differ as markedly as the uncer- 
tainties we assume. Likewise, the dust-to-gas ratio will 
only depend on the condensation fraction in the wind, 
which should not vary as much as a factor of 5 x between 
stars at the same metallicity. This reduces the overall 
error among stars in the same cluster to a factor of ~2x 
and allows us to make direct comparisons between stars. 

3. DISCUSSION 

The wind parameters listed in Table [T] pose several 
problems for stellar wind theory, which we discuss in this 
Section. Firstly, the projected terminal velocities we de- 
rive are very low, implying radiation pressure from the 
star is not effective in driving the wind: we must consider 
whether other processes may contribute to, or dominate, 
the wind acceleration. Secondly, when one considers the 
evolutionary timescales of the stars in question, our de- 
rived mass-loss rates imply that the stars are losing much 
more material than is predicted from stellar theory. Fi- 
nally, we can use these problems to constrain some of the 
wind parameters listed in Table [51 

Mass-loss rate, pulsation amplitude and pulsation pe- 
riod are all known to correlate with stellar evolution. 
This is demonstrated in Figures [SJ [3| and [4] These can 
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Temperature (K) 

Fig. 2. — Hertzsprung-Russell diagram, showing only the upper- 
most parts of the RGB and AGB. Known long-period variables are 
circled, with circle sizes corresponding to their mass-loss rates. 

be used to help solve the above problems. 

3.1. Efficiency of wind acceleration 

A major problem when calculating mass-loss rates from 
metal-poor stars is that the terminal velocities predicted 
from radiative transfer codes are very low. This is a par- 
ticular worry with DUSTY, as no correction is input for 
stellar gravity. In this situation, it is impossible not to 
have mass lost from a star when dust is present, whereas 
the real situation may be that radiation pressure on dust 
is insufficient for it to overcome the stellar gravity field. 
Assuming the dust layers around these stars do have 
some outward velocity, a different method of wind ac- 
celeration may be required. 

Dusty provides a minimum mass for a radiative trans- 
fer model. If a star is less massive than this, the er- 
ror due to the unspecified correction for stellar gravity 
will be greater than the 30% inherent in the model. As 
this error increases to 100%, it is unclear whether radi- 
ation pressure is capable of driving the wind. This is a 
particular concern for the less-luminous stars with high 
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Fig. 3. — Top panel: pcfiod-AIn's diagram for 47 Tuc's 
long-period variables, overplotted with relations derived for the 
Large Magellanic Cloud ( ,Ita et al. 2004): bottom panel: period- 
luminosity diagram for the same stars. Dot size corresponds to 
mass-loss rate. 

mass-loss rates: x03, V18 and V13. While expansion 
velocities of 3-4 km s~^ are not unreason able when com - 
pared to Halo carb on stars (|Groenewegen et al.l Il997i : 
iLagadec et al.|[2010l ). the lower velocities we derive for 
those stars losing less m ass become increa s ingly unre- 
liable due to this error. IBowen fc WilliH^ (fTOOl pre- 
dict that stars with metallicities below [Fe/H] = -1 can- 
not drive winds via radiation pressure, and must rely on 
acoustic (pulsation) d riving to eject mass from the star 
(see also lBowen|[l988[ ). At [Fe/H] = -0.7, 47 Tuc is near 
this hypothetical limit, meaning we must consider accel- 
eration by pulsation and other sources of energy. 

The velocity amplitude of the photosphere in these 
variables ranges from 4-23 km s~^, with VI, V2 
and V3 showing line doubling in near-IR and/or line 
emission in the optical due to the propagation of 
shocks in the atmosphere ( Table [TJ iLebzelter et alj 
[20051: [McDonald fc van LoonI [2001 . While this veloc- 
ity greatly exceeds our calculated wind velocity of ~4 
km s~^, the bulk material does not necessarily attain 
this velocity. Nevertheless, the presence of shocks in 
the extended atmosphere implies that energy is being 
transferred to larger radii, and implies that this energy 
is available for driving mass loss. 

A comparison between the presence of circumstellar 
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Fig. 4. — Mass-loss rate versus a variety of stellar p aramet ers. 
Error bars of a factor of two are added, as discussed in il2.1.5l 
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dust and pulsation reveals some interesting results. All 
the stars listed in Table [T] are variable: Lee 1424 — 
ACVS 002331-7222.6, while x03 shows a variation of 
bKg « 0.15 mag (N. Matsunaga, private communica- 
tion). Conversely, 22 of 43 variable stars show at least 
probable evidence of IR excess. Dust-producing stars 
therefore appear to be (in 47 Tuc) exclusively variable, 
and variable stars tend to be dust producing. In both 
Figures [2] & [3] there appears no easy way of discerning 
RGB-pulsators from AGB-pulsators: the outliers in Fig- 
ure [2] are such due to crowding in the cluster core (Pa- 
per III), while the pulsation periods in Figure [3] do not 
show the bimodality that one might expect from hav- 
ing RGB and AGB stars at different masses. However, 
a bimodality may exist in mass-loss rate with a tran- 
sition occurring around 2500-3000 L0 (Figure S] cf. a 
similar possible bimoda l ity in the gas-production rate in 
iMcDonald h van LoonI (|2007[ )). The five stars with el- 
evated mass-loss rates are also those with the highest 
pulsation velocity amplitudes and among those with the 
longest pulsation periods. As the presence of pulsation 
and its amplitude both appear to correlate with mass-loss 
rate, an obvious conclusion would be that pulsation en- 
ergy is dissipated in the extended atmosphere and helps 
drive the wind, though we caution that the strongest 
variables are also the most-luminous stars with the high- 
est capacity for driving a wind through radiation pres- 
sure. 

Further evidence that pulsation sets the mass-loss rate 
can be drawn from the apparent invariance of mass- 
loss over the narrow range of temperatures we probe 
(Figure S]). The photospheric temperature sets the 
radius at which dust can condense and therefore (if 
other parameters are invariant) the density of the wind. 
Higher-density winds lead to more-efficient dust forma- 
tion, as the timescale for collision between dust-forming 
molecules is shorter. An invariance of mass-loss rate with 
temperature thus means that the formation of dust has 
little impact on the mass- loss rates present, providing 
further evidence that the mass-loss rate is 'hard-coded' 
into the wind before the dust-formation radius. If this is 
true, the role of radiation pressure on dust may simply 
be to affect the outflow velocity of the wind. 

If the wind is acoustically driven by pulsations, then 
the wind velocity and thus mass-loss rate is liable to be 
several times greater than the values listed in Table [U 
since these sources of energy are not included in our 
model. This poses potential problems, as mass-loss rates 
this high would lead to the dissipation of the star's atmo- 
sphere on much shorter timescales than is inferred from 
star counts. This problem may be solved if one invokes 
different parameters for dust in the wind, which could 
lead to mass-loss rates much smaller than the listed val- 
ues (see §3.5|) . 

3.2. Evolutionary problems with our mass-loss rates 

In Paper HI, we discussed the integrated RGB mass 
loss for a typical cluster star. We noted it was likely 
<0.22 M0 and possibly, to first order, invariant with 
metallicity. This invariance would further suggest that 
the primary driver of RGB mass loss is not related to 
dust production. We had also estimated in Paper HI, 
based on the cluster's luminosity function, that one star 
reaches the AGB tip every 80 kyr. 



These figures, however, do not correlate with our de- 
rived mass-loss rates. We find here that the cluster's ~22 
dust-producing stars are losing a total of R:!l.2 xl0~^ 
Mq yr^^ of dust-containing wind. Above the RGB tip 
(r::200G Lq), where dust production is ubiquitous, 13 
stars are producing 1 x 10~^ Mq yr~^. This implies 
that a staggering RiO.8 Mq is lost in the ~1 Myr which 
a star spends on the AGB above the RGB tip. Clearly, 
this is impossible, as this is the entire mass of the star. 
This argues for a mass-loss rate of order 4-8 times lower 
than we measure. If additional wind driving mechanisms 
are important, this factor will increase further. 

Mass loss via dust-containing winds between 1000 and 
2000 Lq can be similarly determined. Here, nine out of 
the 56 stars in this luminosity range are producing ~2 
xl0~^ Mq yr"-'^, implying 0.16 Mq is lost over a 4.5 
Myr timescale. If one assumes that mass loss from these 
stars is also lower than measured, and that some or all 
of these stars are AGB stars, this implies that mass loss 
via non-dusty winds may be of similar or even greater 
importance in removing the estimated ^0.22 Mq from 
the star between the middle of the RGB and middle of 
the AGB. 

At the rates implied in Table [l] the iron in the at- 
mosphere of VI would be entirely ejected within 10 000- 
40 000 years, depending on the remaining envelope mass 
(0.08-0.3 Mq). V8, V2 and V3 would similarly dis- 
perse the iron in their atmospheres within timescales of 
<90000 years. The Padova and Dartmouth isochrones 
used in Paper HI predict that these stars (especially V2 
and V3) still have several hundred kyr left on the AGB, 
implying that our iron mass-loss rates are (again) a factor 
of several too high. 

3.3. Constraining the wind parameters 



In W6.ll we find radiation pressure may not be the 
dominant driver of winds from metal-poor stars, and is 
perhaps incapable of driving them in some cases. The 
simplest solution is to provide momentum from another 
source, such as stellar pulsations, in order to increase 
the wind velocity. Unfortunately, this also has the effect 
of increasing the mass-loss rate, which is already higher 
than stellar evolutionary theory can accommodate fi j3.2l) . 
While we expect another momentum source to be present 
in the most-metal-poor stars (simply because radiation 
pressure on dust is that inefficient), the limits implied by 
the mass-loss rates imply this is not the case in 47 Tuc. 

A solution to this dilemma should address both these 
problems by decreasing the mass-loss rate without signif- 
icantly decreasing the wind velocity. Examining Table [21 
the simplest change would be to bring the dust-to-gas 
ratio closer to unity, which would decrease the mass-loss 
rate and increase the wind velocity. Alternatively, we 
may consider measures that would increase the wind ve- 
locity without appreciable increasing the mass-loss rate, 
such as changing the size distribution of grains. 

The dust-to-gas ratio, ip, given in Eq. ([3]), ass umes that 
the d ust fraction scales with metallicity (cf. Ivan LoonI 
[2OOOI) . This win hold true if every dust species condenses 
with the same efficiency. Naively, one may imagine that 
the fraction of each element condensed would decrease 
toward low metallicity simply because there is less op- 
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TABLE 3 

Maximum elemental abundances in dust 



Element or 




Relative number density 




compound 




CO 


MgSiOs 


Fe & Ni 


r inai 




star 


forms 


condenses 


condense 


values 


H 


lUUUUUU 


1000000 


1000000 


1000000 


lOUOUOU 


He 


81970 


81970 


81970 


81970 


81970 


O 


155 


59 


26 


26 


21 















U 


Mg 


1/11 

14.1 


14.1 


3.2 


3.2 


U 


ol 


1 1 n 
11. U 


11.0 








U 


re 


6.3 


6.3 


6.3 








INl 


0.6 


0.6 


0.6 








Ti 

_L 1 


U.Uo 


0.03 


0.03 


0.03 


n 
U 




U. 10 


0.15 


0.15 


0.15 


n 
U 


Mn 


0.10 


0.10 


0.10 


0.10 


Q 


Ca 


0.4 


0.4 


0.4 


0.4 


Q 


Al 


1 n 


1.0 


1.0 


1.0 


n 

u 




U 


96 


96 


96 


OK 


lVigt)lU3 


n 
VJ 


11.0 


11.0 


11.0 


i l.U 


AAj3"f 1 1 1 1 vj^Ti 








6.3 


6.3 


6.3 


1 1 1 n 1 /^l/'ol 








0.6 


0.6 


0.6 


Other oxides 














4 2 


Element or 




Relative mass density 




compound 












TJ 
Xl 


1000000 


1000000 


1000000 


1000000 


1 nnnnnn 

iUUUUVJU 


He 


325500 


325500 


325500 


325500 


325500 


CO 





2675 


2675 


2675 


2675 


Others 


5846 


3192 


2100 


1716 


1482 


Total gas 


1331400 1331400 


1330300 


1329900 


1329700 


MgSiOa 








1101 


1101 


1101 


Metallic iron 











352 


352 


Metallic nickel 











35 


35 


Other oxides 














222 


Total dust 








1101 


1488 


1710 












1/778 



portunity for metal-metal reactions to take place in a 
metal-poor environment. From our conventional under- 
standing, it is difficult to see how the dust-to-gas ratio 
in the wind could be substantially closer to parity. 

We can, however, calculate a maximum amount of dust 
per unit gas volume, by taking the major dust-forming 
elements (O, C, Fe, Si, Mg, Ca, Al) plus some other 
potential dust-forming metals (Ni, Ti, Cr, Mn) and de- 
termining how much dust they can form. For 47 Tuc, 
we assume the abundan ces of most elements are scaled 
from the solar values of 'Asolund et al.' f2009) by [Z/H] 
= -0.45 dex, following Alvcs-Brito et al. (2005). We 
take the fol lowing exceptions, however: jFe/H] = -0.7 
dex (lHarri£lll996f ): [He/Hi » -0.02 dex (iDorman et al.) 
[1981 : [0/Fe] = +0.2 dex (|D'Orazi et al.ll2010() : [MeTpb] 
ft! fSi/Fel ^ [T i/Fel = -fO.23 dex and [Ca/Fe] = 0.0 dex 
lAlves-Brito et al. (2005). Table[3]shows the stages of this 
calculation, by which we conclude that tp < 1/778. Here 
we make the simplistic assumptions that all dust-forming 
elements start out in the gas phase, that all carbon forms 
CO, that all silicon forms enstatite (MgSiOs), that iron 
and nickel condense entirely as pure metals, and that the 
leftover metals all form their respective oxides. Chemi- 
cal reactions must be calculated on the basis of number 
density, but these must be converted to mass densities in 
order to calculate a dust-to-gas (mass) ratio. For refer- 
ence, the same calculation using solar metallicities arrives 
atiP< 1/235. 

This approach neglects several points. Carbon may be 



enhanced in these stars due to dredge-ujQ, but unless 
it is sufficiently abundant (C/0 > 0.75) to bind all the 
remaining free oxygen, this does not affect the end re- 
sult. We also neglect water condensation, which could 
use up free oxygen, but this would serve to decrease ^p. 
Forsterite (Mg2Si04) may condense out in preference to 
enstatite, and calcium and aluminium may form into sil- 
icates before enstatite, but these would also decrease tp 
where they affect it at all. The formation of FeO and 
iron- rich silicates would increase ip, but it seems that only 
a small percentage of iron condenses into non-metallic 
forms in these stars. Our stated value of V'max is there- 
fore a relatively robust maximum value. 

3.4. Dust condensation efficiency 

Table [1] suggests that the number density of iron grains 
to other dust species in the winds is very high (typically 
10:1, where it can be measured), but the analysis shown 
in Table [3] would suggest it is of order 1:2 if all dust is 
condensed. This would suggest that the fraction of sili- 
cate dust that is condensed is quite low (~5%), however 
we have not yet accounted for two possibilities: that the 
temperatures of the grain populations may be limited, or 
that the size distributions (specifically the surface-area- 
to- volume ratios) of silicates and iron may be different. 

Differing grain temperatures between dust species may 
exist when the wind is sparse enough that dust-gas col- 
lisions cannot equalise the thermal energy. This would 
mean that more-opaque grains, like iron, are warmer; 
more-transparent grains, like silicates are colder for a 
given radius, and could explain the difficulties we have 
in fitting both the metallic-iron continuum, and the 
10- and 20-/im silicate bands simultaneously f §2.1.4p . 
More significantly, however, if thermal energy cannot 
be equalised via collisions, this will lead to a signifi- 
cant velocity drift between the gas and dust components 
of the wind. This will substantially decrease the pro- 
jected mass-loss rate, while the dust outfiow velocity 
will also increase accordingly. If the gas outflow ve- 
locity is below the turbulent velocity in the wind, this 
may allow gas to preferentially fall back, making the 
star more metal poor. Post-AGB stars and planetary 
nebul ae are known to be metal-p oor for similar reasons 
(e.g. Cami'2002':'Maas et al."2005':'De lgado Inglada et"all 
2009; Johnson & Pilachowski 2010), but it is not clear 
that this process is happening while the stars are still on 
the AGB where the outflow velocities are lower and the 
wind densities at the dust-formation radius are greater. 

Conversely, if dust grain temperatures are equalised 
and the formation temperatures of silicates and iron are 
roughly similar, we may examine the effect of different 
grain sizes by taking two limiting cases. If we take the 
limiting case that both silicates and iron grains follow 
the same continuous distribution of ellipsoids, the sil- 
icate condensation fraction is ~5%, and we arrive at a 
value of t/i ~ 1/2000 to 1 /3000 if all the iron is condensed. 
This would increase our derived Al values by ~60% while 
decreasing Voo by ~30%. This is the opposite effect to 
that needed to reconcile our observations with evolution- 
ary theory, therefore we consider this scenario unlikely. 

^ While efficient dredge-up is normally considered to be limited 
to more-massive stars (> 2 Mp)), carbon-enhanced sta rs of ~0.8 
M0 exist in globular clusters (e.g. Ivan Loon et al]|2007l 'l. 
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In the converse limit, silicate and iron dust has con- 
densed efficiently, but the iron grains have average radii 
20 X smaller. This increases M by a moderate (30-60%) 
amount, but increases foo by much more (100-200%), de- 
pending on the exact size distribution. If one is prepared 
to accept that iron grains grow preferentially in one di- 
rection, e.g. if magnetic effects are important, need le-like 
cylinders would be created fcf. lKemper et al.ll2002D . This 
would lead to a larger Voo, but in this case without an 
increased M. However, if identical fractions of silicates 
and iron are condensed, the length/radius ratio of the 
iron grains would need to be a staggering ^1200, mean- 
ing that the grains must be smaller as well. For example, 
a population of grains 5x smaller with a length/radius 
ratio of 10 allows a similar fraction of silicates and iron 
to be condensed with only a slight increase in M. 

There is no obvious mechanism that would cause 
iron to be preferentially condensed at low metallicities, 
though magnetic effects and/or the basic problem of 
forming grains may be important (an iron grain only 
needs two iron atoms to bond, whereas a silicate grain 
requires a series of chemical reactions to form MgSiOa). 
Higher partial pressures of iron may also assist here, 
though it is not clear why this would happen in a metal- 
poor environment. 

Note that in all cases above, we have assumed grains 
are homogeneous: our calculation of the wind parame- 
ters in i j2.1.4l assumes all grains are a uniform mix of 
the included species, while our discussion above takes 
the opposite limit: where grains of each species exist 
separately in pure form. In reality, grains aggregate 
and form inclusions inside each other, which complicates 
the above arguments. In general, however, the evolu- 
tionary constraints provided in ^13. 21 would suggest that 
dust condensation efficiency is high compared to solar- 
metallicity stars (suggesting ip > 1/1000), and the iron 
grains are smaller and/or more-elongated than both the 
silicate grains and the assumed MRN distribution. A 
combination of these effects is likely required to reduce 
M by the factor of 4-8 required without reducing the 
outflow velocities, Uoo, we derive (Table [T]). 

3.5. Composition and radial driving of dust 

Where we are able to determine it, the composition 
of circumstellar dust and its variation with lumin osity is 
similar to that we see in other clusters in McDonald et al.l 
(|2010( ) (see also Figure [5|). The evolution of dust compo- 
sition in 47 Tuc broadly agrees with that seen in other 
clusters at this metallicity. Metallic iron appears to 
be ubiquitous around dust-producing stars, while amor- 
phous silicates are produced by most stars without any 
obvious dependance on luminosity. 

The published spectra o f several stars (¥4, V8 , V13, 
and weakly in V21; see iLebzelter et al.l |2006[) show 
sharp features a r ound 9.7, 11.3, 13.1, 20 and 32 ^m. 
iMcDonald et al.l ()2010[ ) attribute the 20-/xm feature to 
iron-rich wiistite (Mgi-^Fe^^O), where x ~ 1). Recent 
work by Niyogi, Speck & Onaka (submitted) suggests 
that all these features, with the exception of the 13.1 
fxm fe ature (which is pro bably due to an aluminium ox- 
ide — ISloan et all[200l . can be explained by iron-rich 
crystalline silicates. Here, the 20-fim feature is still pro- 
duced by an Fe-0 bond, but it exists in a large crystalline 
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Fig. 5. — Dust type and mass-loss rate as a function of luminosity, 
showing additional data sourced from McDonald ct al. (2010) and 
Paper I. Filled, red squares have only metallic iron dust, hollow 
blue circles also show silicates, while crosses show the stars with 
obvious 13.1 and 20-^m features. Dot size corresponds to estimated 
mass-loss rate. 

lattice of fayalite (Fe2Si04) or ferrosilite (FeSiOa). 

Iron-rich silicates are not expected to form in cir- 
cumstellar environments, as their condensation tempera- 
tures are below that of both magnesium-rich silicates and 
metallic iron (Gail & Scdlmavr 1999). If iron-bearing 
crystalline silicates are a significant component of some 
winds, they must exist alongside metallic iron and amor- 
phous silicate species. One way this could occur is if 
cooling of the wind proceeds at a faster rate than dust 
condensation, which may occur in the tenuous winds of 
low-luminosity, metal-poor stars. Figure [5] shows that 
globular cluster stars which show a 20-/im feature may 
be confined to luminosities below ^4000 Lq. It is diffi- 
cult to say much about their metallicity dependance due 
to insufficient data, especially given silicate production 
appears to decline below [Fe/H] ~ -1. 

Interestingly, the winds found in such stars (x03, V13 
and VIS) appear to be quite cool. In more-luminous 
stars, we find wind temperatures to be higher (~1000 
K). The condensation of metallic iron before silicates 
remains a mystery: for pressures where iron condenses 
around ^1000 K, magnesium-rich sili cates should have 
already formed closer to the star (e.g. IGail fc SedlmavJ 
[1999, their Fi gure 2). Significant departures from LTE, 
or unusually-high partial pressures of iron (caused, per- 
haps, by magnetic confinement or clumping in the wind) 
may offer ways to circumvent this problem, but to invoke 
them as solutions would be speculative without further 
data. 

By whichever process metallic iron is created, if radi- 
ation pressure is effective in wind driving, metallic iron 
condensation should have a dramatic effect on the wind. 
Its condensation should cause a rapid acceleration of the 
wind, dropping the density and preventing other dust 
species from condensing. This may be what we see in 
moderately-luminous stars, such as V3, LWIO and VI 1. 

If radiation pressure is not an effective driver, the gen- 
eral absence of other dust species may be explained if sil- 
icates form using iron grains as seed nuclei. In this case, 
denser or more-metal-rich winds effectively condense sil- 
icates onto iron grains via collisions, while less-dense or 
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metal-poor winds expand before much condensation oc- 
curs, leaving bare iron grains. This scenario may explain 
the iron grains found in the interplanetary dust parti- 
cles know n as GEMs (glasses with em bedded metal and 
sulphides; [Bradley 1991; lMartinl[T995h . 

In the most-evolved stars (VI, V8, V2 and V4), how- 
ever, we do see silicate emission making an increasingly- 
strong presence. Meanwhile, there is a hint that the 
temperature in the dust formation zone (as traced by 
metallic iron condensation) may be decreasing. As noted 
in §3.41 different dust species may also have different 
temperatures. More detailed modelling of the winds 
of these stars, allowing for thermal decoupling of the 
grain species, is required if we are to speculate on factors 
which cause this temperature change. It could, however, 
be linked to stronger pulsation of the star and opacity 
shielding within the wind. In a cooler wind, however, sil- 
icates once again condense out in increasing preference 
to metallic iron, which contradicts our tentative finding 
that the silicate grains are cooler than the iron grains in 
these stars. 

4. CONCLUSIONS 

In this work, we have used spectral energy distributions 
and infrared spectroscopy to establish mass-loss rates 
and investigate mineralogy around the dust-producing 
stars in 47 Tuc. We summarise our conclusions as fol- 
lows: 

• Mass loss in the cluster is dominated by its most lumi- 
nous stars (VI, V8, V2, V3 and V4), however substan- 
tial mass loss appears to be taking place in moderate- 
luminosity stars too (x03, V18, V13). 

• The total mass-loss rate is estimated to be ^1.2 x 
10~^ Mq yr~^, mainly in the form of metallic iron, 



but this is probably over-estimated by a factor of 4-8 
due to unconstrained properties of the dust grains. 

• It remains unclear whether radiation pressure on dust 
is capable of driving a wind alone, due to the low pro- 
jected terminal velocities of the stellar winds. The ex- 
istence of circumstellar dust appears to correlate with 
the presence of stellar pulsation. We theorise that such 
pulsation may inject energy into the wind acoustically. 

• A dust-to-gas ratio closer to unity and smaller or more 
needle-like iron grains are suggested as methods to de- 
crease the total mass-loss rate while allowing for an 
increase in the terminal velocity. 

• Variations in the wind chemistry b r oadly concur with 
those found in [McDonald et all ()2010[ ). with sili- 
cates found preferentially in high-luminosity, high- 
metallicity winds, while iron is found preferentially in 
low-luminosity, low-metallicity winds. Iron may con- 
dense preferentially at low metallicity. Alternatively, 
iron grains may be coated by silicates in dense, metal- 
rich winds, but this coating process may not be effec- 
tive in winds around lower-luminosity and metallicity 
stars. 
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